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In this work we investigate the betweenness centrality in geographical networks and its relationship 
with network communities. We show that vertices with large betweenness define what we call 
characteristic betweenness paths in both modeled and real-world geographical networks. We define 
a geographical network model that possess a simple topology while still being able to present such 
betweenness paths. Using this model, we show that such paths represent pathways between entry and 
exit points of highly connected regions, or communities, of geographical networks. By defining a new 
network, containing information about community adjacencies in the original network, we describe a 
means to characterize the mesoscale connectivity provided by such characteristic betweenness paths. 
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I. INTRODUCTION 

As a consequence of their potential to represent real 
systems, complex networks have become the subject of 
growing interest in recent years [1]. Some examples of 
systems that have been studied using complex networks 
are electric power grids [5], airline routes [3], the World- 
Wide Web U, among others [S]. Networks whose vertices 
are embedded in a space with well-defined coordinates 
and their connectivity is related to the distance between 
the vertices are named geographical networks j^. Ex¬ 
amples of geographical networks include street, airline, 
and neural networks. Models to generate artificial net¬ 
works with geographic dependency have been created |7] 
in order to better understand the real systems. 

The spatial influence on the connectivity gives rise to 
many interesting characteristics |8]. In a previous work 
about the influence of gene expression on neuronal char¬ 
acteristics [S], we computed the betweenness centrality 
m of neuronal networks. While visualizing those neu¬ 
ronal networks we observed that the vertices having the 
highest betweenness centrality values were connected to¬ 
gether creating chains of vertices, which here we call be¬ 
tweenness paths. The presence of betweenness paths in 
the neuronal network suggests a two-level hierarchy of the 
system structure, where the betweenness paths form the 
main pathways for information flow and the remaining 
paths distribute this information locally. Such hierarchi¬ 
cal organization has also been observed for road networks 
[TTHTS]. In [TB] its argued that road networks should be 
planned in different scales, according to their purpose, 
function or management policies. Such scales range from 
the entire city level down to local streets used for prop¬ 
erty access or pedestrian movement. According to [TlllIH] 
the traffic in road networks, modeled by applying the be¬ 


tweenness centrality, and the movement of people in a 
city have also been shown to be power-law distributed, 
meaning that a few paths concentrate most of the city 
traffic. 

Given that a hierarchical structure of pathways has 
been observed in a range of systems, it is interesting to 
study the reasons why they are formed, their importance 
for the correct operation of the system, their influence 
in the network flow, and in which kind of geographi¬ 
cal networks the betweenness paths can be found. To 
better understand the hierarchical structures, it is nec¬ 
essary to study how the network topology is inducing 
those paths. Furthermore, the knowledge about between¬ 
ness paths and their characteristics can be used for better 
planning the system. 

In this work, we study street networks composed by 
two hierarchical levels: a main structure forming the 
backbone of the network, and secondary connections sup¬ 
ported by such backbone. We considered networks gen¬ 
erated from geographic models and real-world street net¬ 
works. From these networks, we identified connected 
vertices with high betweenness centrality. We observed 
that these vertices provide a good covering of the topol¬ 
ogy, in the sense that any network vertex is aways near 
at least one or more characteristic betweenness paths. 
We found that such paths tend to be related with the 
community structure m of the network. So, to bet¬ 
ter understand the characteristic betweenness paths and 
its relationship with communities, we created a random 
geographic model and used community detection to di¬ 
vide the space in well-defined regions. In order to study 
the relationship between such communities and their ge¬ 
ographical positions, we created a sub-summed commu¬ 
nity network, where each vertex is a different community 
and the connection weights are based on the amount of 
shortest paths crossing the community border. 
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This paper is organized as follows. Section |TT] presents 
the basic network theory concepts that are used in this 
paper. In Section |III| we discuss betweenness paths in 
network models and real cases. In Section EYl we create 
a network model to better understand the betweenness 
paths. In Section|V] we describe the relationship between 
geographic and topological measurements. 


II. BASIC CONCEPTS 

There are many ways to compute network centralities, 
each being developed to evaluate different network char¬ 
acteristics [7]. 

In order to locally quantify the network centrality, it 
is possible to measure the degree ki for each vertex i. 
Furthermore, in weighted networks, we can consider the 
edges to define a measurement called strength j8] , which 
is computed for each vertex i as 

= ( 1 ) 

jeri 

where is the set of neighbors of vertex i and Wij is the 
weight of the edge connecting i to j. 

The network centrality can also be computed globally, 
that is, considering the whole network structure. One 
important global centrality measure is the stress central¬ 
ity (SCu) [2^. This measure is defined for each vertex u 
as 


SCu = Y^(7{i,u,j), ( 2 ) 

ij 


where a{i,u,j) represents the number of shortest paths 
between the vertices i and j crossing the vertex u. If 
one vertex is removed from the network it can affect the 
stress of all other vertices. 

Another important global measurement is the between¬ 
ness centrality, which also considers the amount of short¬ 
est paths crossing a vertices [10] . This measure is defined 
as 
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where a{i,u,j) represents the number of shortest paths 
between the vertices i and j crossing the vertex u and 
is the total amount of shortest paths between i 
and j. The betweenness centrality can be implemented 
using a fast algorithm Eij na, where the betweenness 
centrality is computed and the result is used to compute 
the stress centrality. 

Betweenness centrality is frequently used to study 
flows in networks [231 - 125] and can be used in several appli¬ 
cations. In social networks it can be described in terms of 
its message communications, considering that each mes¬ 
sage take the shortest path [7^. In power grid networks. 


this measure can quantify the load of each vertex [Hj. It 
is also used in the cascading-failure dynamics |28| . 

In addition to the centrality measures, another impor¬ 
tant concept studied in complex networks theory is the 
community structure of the system. Where a community 
is defined as a set of vertices densely interconnected, and 
with few connections to vertices outside this set. This 
property can be found in many real networks such as so¬ 
cial networks [29], metabolic networks [30], and in the 
structure of cities m- In city networks, those commu¬ 
nities are shown to be subdivisions of urban structure 
emerging urban elements and their interactions m- 

A common approach to quantify how well-defined are 
the communities in the network is the modularity mea¬ 
sure |19| . The modularity is represented as 

<3 = ^ E - ^) (4) 

where m is the number of network edges, A is the network 
adjacency matrix, ki and kj are the vertex degree i and 
J, respectively, ci and Cj are the communities of vertex 
i and j, respectively, and S{ci,Cj) is 0 if Ci ^ Cj or 1 if 

Ci — Cj. 

Network communities can be detected through the op¬ 
timization of modularity. For a comprehensive review of 
modularity optimization methods, please refere to |32| . 
In this work we use the fast greedy [33] algorithm to de¬ 
tect the communities. 


III. CHARACTERISTIC BETWEENNESS 
PATHS 

In this section we define and analyze the character¬ 
istic betweenness paths. We compute the betweenness 
centrality of all the vertices. Most of these vertices are 
connected, creating chains of vertices. Thus, most of 
them have degree two. We define a betweenness path as 
a sub-network considering only the the highest between¬ 
ness centrality values. In other words we selected only 
vertices with betweenness centrality higher or equal than 
a threshold and created a sub-network composed by such 
vertices. In order to do so, we compute the betweenness 
centrality of the vertices in a network and divide them 
into two groups. Vertices having betweenness larger than 
a threshold T are called characteristic vertices, while the 
rest of vertices in the network are called secondary ver¬ 
tices. 

We will show four networks in which the characteristic 
paths can be observed. Two of these networks are artifi¬ 
cially generated, and the other two are street networks. 
In the first model, known as Waxman [34], the vertices 
are organized randomly in a plane. For each vertex, the 
probability to connect with all other vertices follows an 
exponential decay of the distance between them. In the 
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second model, known as random geometric graph |35| . 
the vertices are randomly organized and the connections 
are created only when the geographical distance is lower 
than a fixed value. We compared the network models 
with two real street networks: from San Joaquin County 
in USA and Oldenburg city in Germany |36| . 

In order to find a reasonable value for the threshold, 
we computed the fraction of characteristic vertices ac¬ 
cording to some selected values of betweenness threshold. 
The fraction is found by dividing the number of charac¬ 
teristic vertices by the size of the network. We perform 
these tests using the street networks, the results can be 
seen in Figure [l] Moreover, we also plot in Figure [l] 
the largest connected component (giant component) of 
characteristic vertices. The majority of the characteris¬ 
tic vertices are included in the giant component for all 
different threshold values. Furthermore, the curves are 
similar for both cities. There is no indication of a proper 
threshold for defining the characteristic vertices. There¬ 
fore, we empirically selected a threshold for which 20% 
of the network vertices become characteristic vertices. If 
we consider the subnetwork formed by such vertices and 
their respective connections it is possible to show that 
for each analyzed network an interesting structure is ob¬ 
tained. Such structure contains a chain-like connectivity 
between the characteristic vertices, meaning that the ma¬ 
jority of vertices have degree two. Such chain-like struc¬ 
ture contains most of the flow in the network, represented 
by the betweenness, and can be seen as the backbone of 
the original network. The secondary vertices rely upon 
this backbone for proper communication with the rest of 
the system. The characteristic betweenness paths were 
also found in other works, but their origins were not an¬ 
alyzed miiisiiiT]. 

In order to illustrate the betweenness paths, we simu¬ 
lated the Waxman network with 1707 vertices and aver¬ 
age degree 2.20 (The network is shown in Figure]^ a)). 
The random geometric graph has 4883 vertices and av¬ 
erage degree 6.18 (as can be shown in Figure |^b)). The 
street network from San Joaquin has 14503 vertices and 
average degree 2.75, and the network representing Old¬ 
enburg streets has 2873 vertices and average degree 2.57. 
The Oldenburg network is shown Figure|^c), where each 
color represents a different community. The San Joaquin 
network is not shown in Figure because the network 
is very large, which hinders a proper visualization of its 
paths and communities. 

For each network we computed the betweenness cen¬ 
trality and highlighted the 20% of vertices with largest 
betweenness values, which defines the characteristic be¬ 
tweenness paths. The characteristic betweenness paths 
occurred in all examples, which are the characteristic ver¬ 
tices. By visual inspection, in each network the charac¬ 
teristic vertices seem to define a backbone structure hav¬ 
ing many two-degree vertices. Such backbone is formed 
by what we call the characteristic betweenness paths. In 



FIG. 1. (Color online) The plot shows the fraction of charac¬ 
teristic vertices against the betweenness threshold. For each 
network we also show the number of characteristic vertices in 
the giant component. 

street networks this result can indicate a tendency that 
the volume of cars is bigger in some streets than in oth¬ 
ers if we consider that all drivers are using the shortest 
paths for the formation of high flow streets. 

In every example we could observe the characteristic 
betweenness paths and it was possible to identify more 
than one component with high betweenness centrality in 
the same network. The component sizes were different for 
each network. In the random geometric model, the giant 
component represents 85% of all vertices in the charac¬ 
teristic path. In the other networks the giant component 
represents more than 98% of all vertices in the paths. 

We analyze the coverage of the characteristic paths. 
In order to do so, we measured the shortest distance 
between each secondary vertex i and preferential paths 
(represented as di), because the shorter the distances 
are, the better are the paths coverage. In order to pro¬ 
vide a fair comparison between different networks, the 
calculated distances are normalized by the average path 
length, I, of each network. Therefore, we define the nor¬ 
malized shortest path distance, 0i, as 



We generate 60 samples of each model. In this case, we 
fit the model variables in order to compare with the Old¬ 
enburg network (modularity 0.92). The samples of Wax- 
man model have 2976.70 ± 16.61 vertices, average degree 
2.71 ± 0.04, and modularity 0.79 ± 0.01. The samples of 
random geographic model have 2892.33 ± 58.46 vertices, 
average degree 4.50 ± 0.09, and modularity 0.97 ± 0.001. 
In order to do a comparison considering the statistical 
variations of the models, these results are shown as his¬ 
tograms in Figure We observed that the fraction of 
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(c) 

FIG. 2. (Color online) Geographical networks, where black circles and highlighted edges show the preferential paths. Item 
(a) represents the network of Waxman model, item (b) represents the network of the random geometric model, and item (c) 
represents the network of Oldenburg streets, where each color represents a different community. 
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vertices with the fewer distance (first bin in the his¬ 
togram) is longer in the random geometric model than 
the Waxman model comparing the two models, we ob¬ 
serve that secondary vertices tend to be closer to charac¬ 
teristic paths in the random geometric graph model, but 
in this network the longest distance is small than that 
case (as can be seen in Figures a) and (b)). Further¬ 
more, the standard deviation is smaller in two networks 
for both models, showing that this effect is not a statisti¬ 
cal artifact caused by fluctuation. The shortest distances 
were also analyzed for the street networks (as can be seen 
in Figures [^c) and (d)). The dashed lines represent the 
average values of 0, which we call the coverage of the 
characteristic paths. The coverage of Oldenburg is larger 
than for San Joaquin, which may be a consequence of the 
higher complexity of Oldenburg’s topology |36| . 

IV. CHARACTERISTIC PATHS MODEL 

In order to better understand the emergence of charac¬ 
teristic paths in geographic networks, we define a model 
able to generate networks with a simple lattice-like topol¬ 
ogy, but still containing a well-defined backbone struc¬ 
ture. The initial structure of the network is a 2D regu¬ 
lar lattice, where the vertices are organized in lines and 
equally spaced columns. First we defined the size, L, of 
the network (number of lines and columns) and two vari¬ 
ables, the maximum distance to connect two vertices as 
dMAX and a probability x to remove randomly the edges 
for each vertex. The algorithm follows three steps: 1- 
The vertices are organized as a square grid with toroidal 
boundary condition; 2- For each pair of vertices i and j 
when the distance dij is less or equal than dm ax an edge 
connecting i to j is created; 3- Each edge is removed with 
probability x. 

In order to visualize the betweenness paths in the 
model presented in this section, we created a network 
using the parameter L = 200, so that the network 
has N = 40000. The connectivity range was set to 
dMAX = v^- Therefore, each vertex of the network 
has degree 8. We randomly remove the connections with 
probability x = 0.65. 

We measured the betweenness centrality for each net¬ 
work vertex and the results can be seen in Figure]^ a). 
We can observe vertices with higher values (characteris¬ 
tic vertices) generating preferential paths in the network, 
in a similar manner as main streets. Our hypothesis is 
that the network communities are inducing the forma¬ 
tion of betweenness paths through connections between 
them. In Figure |^b) we show the network communities 
and characteristic paths (indicated in black). The result 
suggests that the paths are going through communities. 

To compare the characteristic betweenness paths dis¬ 
played by our network model and real networks, we exe¬ 
cuted the same tests as in the previous section. In order 



(a) 



(b) 



(c) 



(d) 


FIG. 3. The histograms show the normalized shortest path 
distance, 0, between each secondary vertex and preferential 
paths. We consider the (a) random geometric model, (b) 
Waxman model, (c) the Oldenburg network and (d) the San 
Joaquin network. The dashed lines represent the coverage and 
error bars in items (a) and (b) shows the standard deviation. 


to infer the statistical variation of the model, we created 
60 network samples. We used the parameter L — 56 
and we randomly removed the connections with proba¬ 
bility X = 0.7. These parameters were set so that the 
generated networks have a size and average degree as 
close as possible to compare with the Oldenburg net¬ 
work. The network measurements are: number of ver¬ 
tices 2829.95 ± 21.34, average degree 2.59 ± 0.03, and 
modularity 0.94 ± 0.002. 
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FIG. 4. (Color online) Item (a) presents a heat map showing the betweenness centrality of vertices belonging to a simple 
realization of our model. The color blue shows regions without vertices. In item (b) each color represents a different community, 
some distinct, non-adjacent, communities are associated to the same color due to color palette limitations.. The color black 
shows the preferential paths. The color white shows regions without vertices. 



FIG. 5. Histogram of normalized shortest path distance, 0, 
between each secondary vertex and preferential betweenness 
paths. 



FIG. 6. Histogram of the coverage for 60 generated networks 
from each considered model. The dashed line represents the 
coverage in Oldenburg’s network. 


Our model has similarities and differences from the 
street networks and the other models. The giant com- 
ponet represents (85 ± 8)% of all characteristic vertices. 
Regarding the shortest distance between secondary ver¬ 
tices and preferential paths (Figure]^, the secondary ver¬ 
tices tend to be more distant from the backbone than in 
the Oldenburg and Waxman networks. 

We compared the average of shortest distance between 
the secondary vertices and the preferential betweenness 
paths in network models and the Oldenburg network. For 
each generated network we calculate coverage, which we 
represent by (0), and in Figure we show the distri¬ 
bution of (0) for the generated networks. The dashed 
line shows the coverage for Oldenburg. The coverage for 
Oldenburg is closest to the measured values for the char¬ 
acteristic path model. Besides, the standard deviation is 
low, which confirms that the betweenness paths have a 
similar structure in all generated networks. Therefore, in 
the following we further analyze the betweenness paths 
using a single realization of our model. 


To quantify our network communities in terms of their 
respective degree and betweenness centrality, we con¬ 
sidered two types of vertices, the geographical border 
community vertices and the internal community vertices. 
The first type represents vertices that have at least one 
neighbor in the grid belonging to a different community, 
no matter if they are connected or not. The second type 
represents a vertex that is not a geographical border. In 
order to define each vertex type, we considered all eight 
geographical neighbors, of the vertices. In some analy¬ 
ses we also considered the topological border community 
vertices, which are the vertices connected with one or 
more vertices from different community. We generated a 
network with L = 200 and average degree 2.92 according 
to our model. Using this network we computed three dif¬ 
ferent comparisons. In the first, we compare the average 
degree between internal vertices and border community 
vertices (as can be seen in Figure]^ a)). In this compar¬ 
ison it is possible to observe that the average degree in 
the internal community vertices is higher than the border 
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community vertices. Thereby, we can affirm that the ver¬ 
tex degree is influencing how the community is formed. 

In order to describe how the vertices are connecting 
communities, we measured the average degree in the ge¬ 
ographical community border and in the topological com¬ 
munity border, represented in Figure |^b). In 98% of the 
communities, the average degree of vertices belonging to 
the geographical border is lower than for vertices con¬ 
tained in the topological border. 

We made a similar comparison, between geographic 
and topological borders, but now measuring the between¬ 
ness centrality instead of the degree (as can be seen in 
Figure [^c)). This analysis can show if connections be¬ 
tween communities (topological border) are generating a 
characteristic betweenness paths. Figure shows that 
the average geographical border betweenness is lower 
than the topological border betweenness for all commu¬ 
nities. Thus, in average, the vertices where the paths 
are crossing the community borders have higher between¬ 
ness centrality and degree than the geographical borders. 
Some vertices of characteristic paths are part of the topo¬ 
logical border. In other words, they are connecting the 
communities. 

The results indicate that the paths represent routes of 
efficient communication between the communities. This 
idea is illustrated in Figure Vertices included in the 
highlighted gray region belong to the same community. 
The community has three vertices on its topological bor¬ 
der and the shortest paths connecting them, shown in 
red, are creating the preferential paths with higher be¬ 
tweenness centrality. 


V. COMMUNITY ADJACENCY NETWORK 


In order to observe the relationship between communi¬ 
ties and their importance, considering geographical char¬ 
acteristics, we created a network based on the topolog¬ 
ical adjacency between communities of our characteris¬ 
tic paths model. Through this model it is possible to 
compare the geographical characteristics (e.g., the com¬ 
munity area) with topological characteristics, such as the 
degree of a vertex, which represents the number of com¬ 
munities connected with a given community. In this net¬ 
work, each vertex represents a community. A weight 
equal to SCi -I- SCj (defined in Equation 2) is associ¬ 
ated to each edge connecting communities i and j. The 


resulting network is undirected. In Figure 11 we show a 
typical realization of this network. It is possible to ob¬ 
serve small communities having strong connections, in¬ 
dicating that betweenness paths are also present in such 
small communities. On the other hand, there are big 
communities with weak connections, showing that few 
preferential paths are crossing them. 

We compared local topology characteristics of the com¬ 
munity adjacency network with geometrical characteris¬ 


tics (the scatter plots can be seen in Figure 121. Fig¬ 
ures [I^a), (b) and (c) show the degree as a function 
of area, perimeter and diameter, respectively and Fig¬ 
ures Eld), (e) and (f) show the same measures using 
the vertex strength instead of the degree. The degree 
is highly correlated with the geographical characteristics 
of communities. On the other hand, the strength has 
little relationship with the communities geometry. This 
happens because the strength is uncorrelated with the 
geographical characteristics of communities. 


In order to validate these results, we created an artifi¬ 
cial geographic network with fixed communities. In this 
network, each community is a regular lattice, where the 
vertices are connected with their eight nearest neighbors, 
and the communities are interconnected through some se¬ 
lected vertices. Figure]^ shows the network and highlight 
the highest values of betweenness centrality (the charac¬ 
teristic vertices). In this test, it is possible to observe 
that the characteristic paths are connecting the commu¬ 
nities. Besides, through this model, it is possible to note 
that the paths are not created at random, that is, their 
emergence have a well-defined mechanism. 

It is also interesting to analyze if the number of short¬ 
est paths connecting pairs of communities are related to 
the sizes of the communities geographical border. In or¬ 
der to do so, we plot the perimeter of adjacency between 
communities against the sum of the vertex stress connect¬ 
ing such communities, which can be seen in Figure |10| 
Besides, we calculated the Pearson correlation coefficient 
between the measurements, obtaining a value of 0.23. 
Thus, these measurements are uncorrelated. 


VI. CONCLUSION 

The structure of a geographical network has a direct 
impact on how the system can be navigated. Our analysis 
revealed that betweenness paths are a common charac¬ 
teristic of geographical networks, since they can be found 
in geographical models as well as real world networks. 
We characterized the betweenness paths in terms of the 
distance between each network secondary vertex and the 
closest betweenness path. Such measurement can be re¬ 
lated to the covering provided by the paths to reach the 
entire network. We showed that in both the geographi¬ 
cal models and the street networks the betweenness paths 
provide a good covering of the system. In order to better 
understand the establishment of the characteristic paths, 
we proposed a model to generate geographical networks 
having a simple topology but that can still originate such 
paths. 

It was found that the characteristic betweenness paths 
have a strong relationship with the presence of well- 
defined communities in the networks. We compared 
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FIG. 7. (Color online) Image (a) shows the degree comparison between internal vertices and border community vertices, image 
(b) shows the degree comparison between topological and geographical border, and image (c) shows the betweenness centrality 
comparison between topological and geographical border. The community indexes were ordered according to the degree of 
internal community vertices. 
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FIG. 8. (Color online) The highlighted region represents a 
community and the red squares represent characteristic ver¬ 
tices crossing the community. 

the vertices degrees and betweenness centrality among 
geographical and topological borders of the communi¬ 
ties, as well as the internal vertices of the communities. 
This analysis revealed that the betweenness characteris¬ 
tic paths are usually induced by the communities, in the 
sense that such paths provide a communication between 
entry and exit points of the communities. 

The number of shortest paths passing through the 
topological border of adjacent communities seems to be 
independent of the geometry of the interface between the 
communities. This was revealed by the absence of corre¬ 
lation between the total stress of vertices connecting ad¬ 
jacent communities and the perimeter of the adjacency. 
In order to better understand the relationship between 
communities and characteristic betweenness paths, we 
defined a new structure based on the vertices between¬ 
ness and the community organization of the network be- 



FIG. 9. (Golor online) Network created by connecting lattice¬ 
like communities. Each community is represented by a differ¬ 
ent color and the characteristic paths are highlighted. 
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FIG. 10. Scatter plot showing for each community the sum¬ 
mation of all stress connecting vertices from a community to 
another against the perimeter of adjacency. The Pearson cor¬ 
relation (p) is shown inside the plot. 
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FIG. 11. (Color online) Heat map representing the commu¬ 
nity adjacency network from communities, where the color of 
edges represent the weight (black is zero and white is one). 

ing studied, which we called community network. By 
characterizing this meta structure, we concluded that ge¬ 
ographical characteristics can influence the connectivity 
between communities, which takes place along between¬ 
ness characteristic paths. 

It remains an open question if other, non-geographical, 
network models can also present such characteristic be¬ 
tweenness paths. Also, the proper impact that such paths 
may have on a dynamics taking place on the system is 
also an interesting aspect for future studies. 
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